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ABSTRACT 

The topic of this paper is the stress analysis of an 
icebreaker bow. The purpose in writing this paper was to 
familiarize myself with the various methods that have been 
developed in determining bow loads in ice breaking and the 
methods that have been developed to compute the stresses in 
a ship's structure. From a study of many methods, I wanted 
to select one that would be most useful to me in my day to 
day work when I depart from M.I.T. I selected as the ship 
to model the United States Coast Guard Icebreaker WESTWIND. 

In the first section, the various methods of 
determining bow loads in ice breaking was studied. The 
various methods that were put forth by scholars and engineers 
over the past half-century gave bow loads ranging from one 
hundred and fifty-four tons to values as high as twenty-nine 
hundred and fifty tons of force at the bow in a vertical 
direction. I chose as the best results to use in loading 
the bow, those values which are determinable from the 
computer program developed by Dr. White in the doctoral 
thesis, "The Dynamically Developed Force at the Bo w of an 
;> Icebreaker." The computer program developed by Dr. White 
gives the vertical and horizontal forces acting on the bow 
as function of the length between perpendiculars, the beam 
at the waterline, mean draft, displacement, bow angle, 
spread angle complement, impact velocity, water plane 
coefficient, longitudinal position of the center of floata- 
tion, the longitudinal position of the center of gravity, the 
height of the center of gravity, height of the thrust line 
above the baseline, bollard thrust that would be obtained 
for RPM being used during the crushing phase of the ice, 
longitudinal metercentric height, kenetic friction coefficient 
-*of ice and the compressive failure stress of ice. The verti- 
cal load determined by his computer program for the WESTWIND 
ramming solid ice at fifteen knots gave a peak value of two 
thousand tons . 
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The next section of this paper covers the various 
methods available for determining the stress loading inside 
the structure. The most used computer programs that were 
reported upon were the STRESS and FRAN programs. Both of 
these programs being strictly a frame analysis of the 
structure leave much to be desired in the complete know- 
ledge of the stresses present in a complicated structure 
such as an icebreaker. 

The most promising and most exact analysis that 
is available at the present time, utilizes the finite 
element method. A program developed and presently being 
revised by NASA that quickly and efficiently analyzes almost 
any structure is the "Structural Analysis and Matrix Inter- 
pretive System" (SAMIS) . The remainder of the paper gives 
a short summary of the program, the data required by the 
analyst, and finally the necessary steps required for the 
analysis of the icebreaker WESTVJIND. 
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INTRODUCTION 



In the Shipbuilder by John Ruskin, 1819-1900, he 
states, "Take it all in all a ship of the line is a most 
honorable thing that man as a gregarious animal has ever 
produced. Into that he has put as much of his human patience, 
common sense, forethought, experimental philosophy, self 
control, habits of order and obedience, thoroughly wrought 
handiwork, definance of brute elements, careless courage, 
careful patriotism, and calm expectation of the judgement of 
God as can be put into a space 300 feet long and forty feet 
broad, " (Ref . 1) * 

Since that time, many ships have made history — 
some of them by going to the bottom of the sea; others, in 
the defiance of the brute elements have trespassed the 
Arctic waste. This type of ship is what this thesis is all 
about . 

The icebreaker hull at its birth was merely a 
standard hull with extra reinforcement in the bow. As time 
progressed, changes in design were based largely on the 
success and failure of past designs. Efficient ice breaking 
(or what was believed to be efficient ice breaking) led to 
the thirty degree angle bow. Removal of the unsuccessful 

•k 

References are listed at the end of this report. 
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bow propeller on the wind class breaker produced the step 
which prevented the vessel from riding completely up on the 
ice. With the advent of the computer, more efficient and 
better designs have been produced. 

Although very few papers have been written on the 
subject of icebreaker design, those that have been produced 
have made far reaching effects on the design of the ice- 
breaker hull. 

METHODS OF DETERMINING BOW LOADS 

The earlier works were done by the Russians. 

M. K. Tarshis in his paper, "Ice Loads Acting on Ships," 

(Ref. 2) put forth a formula to determine the impact load on 
the vessel. This formula considers the speed, ,the angle of 
blow, and the square root of the relative mass and the 
relative rigidity of the ship. He uses as an example a 655 
ton displacement vessel in contact with an ice floe that is 
25 meters in diameter and one meter in thickness. From his 
formula he deduces that the impact load is 220 tons. He 
assumed in this work a crushing strength of ice of 570 pounds 
per square inch. 

The general idea of this paper is that the speed 
of the vessel and the angle of blow in the area of contact 
with the ice floe are the major contributors to the impact 
load on the vessel; hence, if you double the speed of the 
vessel you double the impact load on the bow of the vessel. 
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"The paper has a high theoretical content and requires many 
parameters difficult to determine. The relative angles of 
the hull form at the point of impact are necessary and these 
can vary quite appreciably in a relatively short distance." 
(Ref. 2, pg. 5) 

L. M. Nogid in his paper, "Impact of Ships with 
Ice," (Ref. 3) attempts to determine the reduction of speed 
of the icebreaker as it contacts the ice. An interesting 
sidelight of this paper is that the theoretical amount of 
force required to initiate the crack in the ice is quite 
small, but that the force required to propogate that crack 
might be quite large. It also goes on to show that the 
load required to break off pieces of the ice is more than 
three times the load required to initiate a crack. 

In his paper, Mr. Nogid gives a method to determine 
the speed reduction a vessel will have when coming into con- 
tact with an ice floe of given diameter if you know the 
maximum load the ship could withstand at the point of impact 
and the angle formed with the ice at the point of contact. 

The area can also be determined if you know what the strength 
of the ice is in the area in which the ice breaking operation 
is being carried out. 

He divides the forces and the strength requirements 
of the hull into two different areas. One is impact and the 
other is compression. 



3 



Under the impact force, he determines that this 



force is equal to the mass of the ship times the velocity of 
the ship, modified by parameters which are dependent upon the 
angle of the hull at the point of impact, the area of impact 
and other similar items. 

Under compression, the author assumes that the 
loading is directly related to the crushing strength of the 
ice. He also suggests some relationships between the changes 
of length and the strength of the side frames. His formula 
suggests that the strength of side frames is related to the 
cube root of the ratio prototype to the proposed length. He 
also goes further in an attempt to relate hull strength to 
the frame spacing of the vessel. This paper, in general, is 
the theoretical approach to the old idea of taking a success- 
ful design and changing the dimensions to make a more capable 
design . 

U. N. Raskin in his paper, "Method of Determining 
the Stresses in Decks and Transverse Bulkheads Caused by Ice 
Loads," (Ref. 4) uses a compressive ice load on the side 
framing of a ship to determine what the stresses are in the 
decks and bulkheads. He divides the decks into strips with 
the only load carrying area being the deck itself and those 
beams in direct contact with the deck. The stringer is 
considered to be in simple compression, while the remaining 
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strips act as rigid members to bending and carry only shear 
loads. In between these strips are imaginary elastic bands 
which allow relative movement back and forth between the 
bands. The bulkheads are considered as simple plates in 
end compression. 

Using these models, he presents a method to 
determine stresses in these decks and bulkheads. I believe 
it would be rather difficult to accurately model a complicated 
structure, such as an icebreaker, in this form and achieve 
accurate results. 

In September of 1965, R. M. White completed his 
doctoral thesis, "Dynamically Developed Force at the Bow of 
an Icebreaker." (Ref. 5) In this thesis. Dr. White developed 
a computer program for the prediction of the dynamically 
developed force under the bow of an icebreaker while ramming 
solid ice. 

The solution is based primarily on Newton's Law of 
Motion. The problem was broken down into two basic phases. 

The crushing phase represents the local crushing of the ice 
to accommodate the bow. The sliding phase represents the 
sliding up of the bow without further crushing. "The final 
state represents (temporary) equilibrium when motion has 
stopped; the vertical force at the bow at this state is 
relatively sustained and is the most effective in breaking 
the ice." (Ref. 5, pg. 3} 
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Dr. White gives one of the earlier developed force 
equations which is attributed to R. Ruineberg in 1888 (Ref. 6). 



F 



BZ 



T IB^ COS ^B COS ^ sin 

(sin i_ cos 8 + f , cos i„ ) 
B k B 



where 



F n „ = Downward force 
13 Ll 

T g = Thrust available for icebreaking. 

f, = Coefficient of friction (suggested as .05) . 

8 = The angle with respect to the center line plane 

of a normal to the shell. 

i B = The angle the stem makes with the base plane. 



Dr. White states that Ruineberg assumed no momentum 
effects and no forward motion, so that all thrust is applied 
to ice breaking and is applied horizontally at all times . The 
direction of friction force remains the same during forward 
horizontal progress and the trim does not effect the solution. 
For the icebreaker WESTWIND, this would give a downward force 
of four hundred and fifty-four tons for this continuous mode 
of ice breaking. 

The next equation Dr. White discussed was one 
developed by A. Kari in 1921 (Ref. 7) . 
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4480 ACL sin 0 
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where 



A is displacement in tons 
0 is change in trim 

L is the length between perpendiculars 
Ii is the draft in feet 

and 

BM_ x H GM t x H 

__ -Li l_i ^ ~ 



where 

BM^ is assumed to be equal to GM^. 

GM is longitudinal metacentric height in feet. 

-Li 

Kari makes the following assumptions: (1) no 

momentum effect, (2) the bow rises the thickness of the ice, 
(3) the distance from point of contact to the center of 
flotation and the center of gravity are the same, (4) 
effective displacement and the draft is not changed by load 
force, (5) the center of flotation and metacenter remain 
fixed, (6) C = .07, and (7) there is no frictional force. 
This would give a downward force of approximately seven 
hundred and fifty-eight for the icebreaker WESTWIND. 

The third method discussed by Dr. White was one 
done by E. R. Simonson in 1936. (Ref. 8) 

T 

F = in 

BZ tan (iB + 0) 
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In this method, Simonson assumed there are no momentum effects, 
friction is negligible, thrust is directed horizontally, the 
coefficient of friction serves as a pivot point, and there is 
no change in displacement. Dr. White states, "... Simonson's 
equation is limited to being a good approximation for the 
stopped equilibrium position." (Ref. 5, pg. 46) This would 
give a downward force of one hundred and seventy-three tons 
for the icebreaker WESTWIMD . 

The fourth work discussed by Dr. White on the 
downward force developed in ice breaking is one by L. V. 
Vinogradov which was published in 1946. (Ref. 9) Dr. White 
states that this is the first time that force due to ramming 
was put into useable mathematical form. 



o o v 9 

f bz = xt + {x a t a + \r 



v 



[1 - Cl - 



2 • 2. , T7 2 . 
e sm i fi ) - V-J 



* * 



gD 



where 



X = 



1 + 



75 - tan i„ 

cos 3 B 

f k . 

75 - cot i n 

cos g B 



T = thrust in tons 



Y = 



1 + 5 - cot i 

cos 6 B 



Symbols have been changed to previously defined symbols. 
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03 = 


displacement 


D 


draft in ft. 


A = 


^u + kl 1 ] 

C W U k 2 4C W ] 


C B = 


block coefficient 


C W = 


waterline coefficient 


k l = 


q 

L/2 


q 


distance from point of impact to center of flotation 


k 2 = 


A GM L • H • C B 

w 


V o = 


speed prior to impact 


v i = 


speed while sliding up 


change 


Vinogradov assumed that thrust is always horizontal, 
in trim and draft is small enough so as not to affect 



the waterplane characteristics and the metacentric height, 
and GM t = BM_ . 

In Vinogradov's equation the force downward is 
related to the frictional force, the angle of the bow, the 
spread angle of the bow, the coefficient of friction, the 
block coefficient, the waterline coefficient, the distance 
of impact from the center of flotation, the speed prior to 
impact, and the speed while sliding up. This would give a 
downward force of approximately two thousand nine hundred 
and sixty tons for the icebreaker WESTWIND . 
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The fifth equation discussed by Dr. White is one 
by F. Richardson in 1959 in some personal correspondence to 
Dr. White. "The development was almost identical to 
Vinogradov's but did modify some of his weaknesses to some 
extent. For example, Richardson uses a term for the loss of 
energy due to wave and frictional resistance (not ice) from 
the instant of contact up to the moment the ice breaks or 
motion ceases. He also recognizes an effective increase in 
the mass of the icebreaker due to entrained water." (Ref 5,p 49) 
The final equation discussed by Dr. White was one 
developed by V . R. Milano in 1962 which is a modification of 
Vinogradov's equation. Dr. White states, "One of the main 
contributions was to express thrust as a function of 'Bollard 
Puli' ." (Ref. 5, p. 51) 

Dr. White in his analysis assumes: 

1. The force normal to the plating is 
represented by the product of the area 
of contact and the compressive failure 
stress of the ice. 

2. There is a friction force acting in the 
plane of the plating. 

3. The icebreaker is treated as a "solid body." 

4 . The radius of gyration of an icebreaker 
can be assigned 0.266. 
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Crushing has ceased when velocity is in 
the direction of the angle of the bow 
plus the angle of trim. 

6. During sliding phase, the point of 

contact is fixed and is at the level 
of the waterline. 

Results show that the peak load occurs during the crushing 
phase . 

In his computer program, the downward force is 
determined as a function of 



Computer 

Symbol 



L Length between perpendiculars, ft. 

B Beam at waterline, ft. 

II Mean draft, ft. 

A Displacement, lbs. 

i c Bow angle (from base line to stem) , radians 

SA Spread angle complement (normal to bow 

plating with respect to center line 

plane) , radians 
Impace velocity, ft/sec. 

a = C^ a, Waterplane coefficient, dimensionless 



LCF Longitudinal position of the center of 

flotation (- if aft of amidships, + if 
forward), ft. 



BP 

B 

H 

DIS 

BA 

SA 



VI 

AL 

CF 
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Computer 

Symbol 



CG 


LCG, Longitudinal position of the center CG 

of gravity (- if aft of amidships, + if 
forward), ft. 


KG 


Height of center of gravity above base GK 

line, ft. 


d 


Height of thrust line above base line D 

near center of gravity, ft. 


T 

BOL 


Bollard thrust which would be obtained TB 

for rpm used during crushing and 
sliding, lbs. 


GM t 

Li 


Longitudinal metacentric height, ft. GM 


4-1 


Kenetic coefficient of friction of FK 


a 


ice/ship . 

2 

Compressive failure stress of ice, #/ft . SIG 


program 


The following pages are Dr. White's complete 
for developing this downward force. (Ref. 5, 


pp. 385- 


-393) . 
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6 FORMAT (F11.5/3F11.5/3F11.5/3F11.5/2F12.5/4F11.5//) 

IF (XD) 38,38,37 

38 PRINT 39, FZC 
GO TO 36 

39 FORMAT (44H SHIP STOPPED DURING CRUSHING PHASE, FZC2= , E12 . 5//) 
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DISC 4 - WB3-(BE3**2)-(AL3**2)-4.*AL3*AL4-(AL4**2) 
IF (DISC4 ) 23,24,24 

PRINT 13, DISC4 
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FORMAT (22H EXTRACTING THRUST =,E12.5/50H RATIO OF EXTRACTING 
1THRUST TO BOLLARD THRUST IS ,F15.3//) 

GO TO 36 
END 
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STRESS ANALYSIS OF SHIPS STRUCTURES 



In the analysis of the Coast Guard icebreaker 
GLACIER, in Lloyd's Register of Shipping Research and 
Technical Advisory Services Report No. 5095 (Ref. 10), the 
IBM program STRESS was utilized. This is a computer program 
which performs a linear elastic analysis of a framed struc- 
ture under the influence of a static load and is capable of 
solving two or three dimensional problems . 

The condition of loading that was utilized was 
400 psi ice pressure on the midship section of the icebreaker, 
and the ice condition was considered to be 10 feet in thick- 
ness. There were two conditions imposed with this 10 foot 
thickness of ice. 

(1) Two feet of ice extending above the waterline 
and eight feet below it. 

(2) Assume that all the ice was placed at the weakest 
point of the structure that was likely to encounter 
ice load in the ice breaking operation. 

Once again in this program, you have a problem of 
the plating not carrying any load which certainly, in the 
analysis of an icebreaker where you have 2 1/2 and 3 inches 
of plating, will lead to large errors. 

In the analysis that was carried out, two sets of 
end constraints were considered: 
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(1) Supports only at the extreme corners of the 
structure or to have four joints. 

(2) Supports at all inboard ends of keelsons, etc. 
in addition to the condition above going a total 
of eight joints. 

The bending moments obtained from this program with 
four and eight supports varied markedly with factors as low 
as two and with others as high as 60 and none of the varia- 
tions showed any real relationships . The four support condi- 
tions gave values at various positions in the structure that 
at some points were much greater and sometimes much less than 
the corresponding points of the eight support conditions. 

The results of the paper did not give any real 
concrete evidence as to what the actual stresses would be in 
the vessel considered under the 400 psi ice load condition. 
The conclusions were only that the structure should be of 
the grid structure instead of the truss structure. Once 
again, this is nothing more than a selection from the best 
of the successful designs to be applied to a new design. 

Concurrent with the previous report, further 
investigation was also being done on these amidship frame 
areas by Paris Genalis . (Ref. 11) ‘Ie also used the IBM 
computer program STRESS to analyze the amidship 1 s section. 

As he quotes in his report from 
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In his report, he breaks the analysis down into 
four stages: 

(1) Consideration of part of only one frame. 

(2) Consideration of one half of one frame. 

(3) Two dimensional analysis of one complete ring 
frame . 

(4) Three dimensional analysis of more than one 
complete ring frame. 

Each of the analyses carried out was used to 
substantiate the assumptions and methods utilized in the 
next analysis. In the final three dimensional analysis, the 
plating was taken into account by replacing it with stiffeners 
to get away from the problems that Lloyd's of London had in 
their report. 

Each of the problems considered in the three 
dimensional analysis, took a total of one hour or more of 
computer time. The number of calculations was very large, 
hence the probability of error was very large. 

Mr. Genalis also went further in his analysis in 
that he considered machinery weight, steei weight, and the 
effects of varying amounts of buoyancy in addition to the 
ice loads . The wide variations in the results of his analy- 
sis point out the absolute need for proper modeling of the 
structure in order to determine realistic stresses. 



28 



Further limitations of the STRESS program is that 
it cannot handle curved members and thus assumes that each 
member is straight and slender. 

Another application of the IBM program STRESS was 
completed by Lloyd's Register of Shipping, Research and 
Technical Advisory Service Department in Report No. 5051, 

(Ref. 12), in which a structure of shell plating subjected 
to ice pressure was analyzed. "The purpose of this investi- 
gation was to design a simple panel of shell plating stiffened 
in such a way so that the overall and local strength was 
sufficient to withstand a specified ice pressure uniformly 
distributed over the whole panel and for a concentrated load 
acting on one frame only for the full depth of the panel." 
(Ref. 12, p. 1) 

Here again, we have results that do not take into 
account the shell plating other than to simulate it as part 
of the stiffeners. It does give results that compare well 
with simple beam theory, but does not give the information 
necessary for the designers who must, in these days of 
economy, depart from the large factor of safety of five or 
ten to insure that the vessel will be successful and safe 
from the environment it operates in. 

The IBM STRESS program was also utilized by 
Consulter, Inc. of 1725 K Street, N.W., Washington, D. C. 
in "Polar Icebreaker Preliminary Structural Design and 
Special Studies," (Ref. 13) completed in August of 1968, 
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and gave a summary of known ice properties and an analysis 
of possible structures to be utilized in a new type of 
icebreaker called the M-10 at that time. 

The ice loading used was that developed from the 
analysis by M. K. Tarshis, "Ice Loads Acting on Ships," 

(Ref. 14) a translation of a Russian text published in 
Rechnio Transport, Vol . 16, Mo. 12, 1957, pg . 19. 

Using these values of ice loads they determined the 
capabilities of the proposed design along both elastic limits 
and plastic hinge limits. They included factors of 1 1/2 for 
impact loading and two for static condition loading to 
account for the lack of reliability of the values of the 
properties of ice. 

The results once again lack aspects of reality due 
to the limitations of the SJTRESS program. 

During recent years a new method called the finite 
element technique has been developed to give more accurate 
solutions for structure analysis. The first easily under- 
stood, comprehensive, presentation of the method was written 
by 0. C. Zienkiewicz in his boom, "The Finite Element Method 
in Structural and Continuum Mechanics." (Ref. 15) Develop- 
ment of a computer program capable of analyzing a wide 
variety of structures was developed by the Jet Propulsion 
Laboratory, California Institute of Technology, Pasadena, 
California, called the Structural Analysis and Matrix 
Interpretive System (SAMIS). (Ref. 16) 
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THE STRUCTURAL ANALYSIS AND MATRIX INTERPRETIVE SYSTEM 

The Structural Analysis and Matrix Interpretive 
System (Ref. 16) uses elements that are at present restricted 
to flat triangular facets that are joined together along 
their edges, and line elements that are joined to the rest of 
the structure at their ends. All elements are capable of 
resisting stretching, shearing, bending and twisting loads. 
Heating, acceleration and pressure loads can also be analyzed. 
Additional loads can be introduced as energy equivalent con- 
centrated loads at points on the structure. 

Structural changes are defined by stresses, deflec- 
tions, flexibilities and stiffnesses. It also has the capa- 
bilities of computing the natural frequencies of the struc- 
ture and the mode shapes . 

The finite element involved utilizes two ideas 

1. The structure is divided into small elements 

2. The problem is solved by a structural stiffness 
analysis . 

Load deflection relations are taken for each element of the 
structure. The coefficients of these relations form the 
stiffness matrix. When the whole system is taken into account, 
the load deformation relations for the entire system stiff- 
ness matrix is developed by summing the stiffness matrices 
of the pieces composing the system. Where there are common 
grid points, the forces are simply added. Boundary conditions 
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are also formulated in matrix notation. The displacements 
at each of the corners of the facets and the ends of the 



beam that make up the complete structure are determined by 
solving systems of equations simultaneously. By taking 
these displacements, the stresses in each of the elements 
are determined. 

In the Jet Propulsion Laboratory Technical 
Memorandum No. 33-317, (Ref. 17), the mathematical equations 
of the structure are given as: 

[K] {d} - X[K i ]{d} + [C] [d] = {P(t)J - [M s ]{dHa> = [R] { d } 



where {d} = 


The vector of grid point displacements 
(dots indicate time derivatives) . 


(a) = 


The vector of element stresses. 


[K] 


The small deflection stiffness matrix. 


[C] 


The small deflection damping matrix. 


[K.] = 


The initial stress stiffness matrix 


X 


A scalar defining the magnitude of the 
initial stress distribution. 


[MJ = 


The mass loading matrix. 


{P} 


The force loading matrix. 


[R] 


The matrix of stress coefficients. 



Boundary conditions are applied on "d' 1 and "X" to make the 
above equations solvable. 
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The static and dynamic displacement response of 
the structure is given by the first equation and the stress 
associated with a given set of displacements is given by the 
second equation. 

An assumed displacement function which is continuous 
along a facet's edge and continuous over the facet is used to 
develop the stiffness coefficients by minimizing the potential 
energy of the entire structure. The resulting stiffness 
coefficients define the forces and moments at the apexes of 
the triangle which satisfy macroscopic force and moment 
equilibrium conditions. 

The displacements used to develop the stiffness 
coefficients are also used to develop the loading coefficients 
under the restraints of minimum potential energy. 

The stress coefficients are developed using the 
same assumed displacements functions and stress-strain and 
strain deformation relationships. The stress so developed is 
the mean stress at the facet centroid. 

Two coordinate systems are used to implement this 
program. The local system describes the grid points of each 
element as to their local relationship, and the common system 
describes the overall system. In the facet element, the local 
x-y plane is assumed to be coincident with the midplane of the 
facet. The displacements of each grid point are defined as a 
"d" vector which consists of translations in the x, y, and z 
directions, and the rotation about the x and y axis. Each 
component of these deflections is assumed to vary linearly 
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over the midplane surface. The deflection component at any 
point is therefore expressed as 

d . =a.x+b. + c . 

11 ly 1 

In the stress-strain relation, it is assumed that the stress 
in the z direction is zero; assuming that the thickness 
radius ratio is small compared to one for the shell. 

The selection of the triangular elements must be 
such that the stiffness matrix is positive definite. This can 
be done as long as the largest angle for any triangle is less 
than 90 degrees. If an angle is larger than 90 degrees in the 
facet, the stiffness matrix will be indefinite and thus cannot 
be used. The most accurate results are obtained if the 
triangles are all equilateral. 

Discrete loads acting at a joint are described by 
forces acting in the x, y, and z directions and moments acting 
about the x and y axis. When torque is considered, it must 
be treated as a couple with forces in the x and y directions. 

Local coordinate axis are used in developing the 
matrices of coefficient and, hence, the need for a trans- 
formation from the common axis. This is automatically 
accomplished if the grid points are given in the data as 
the common coordinates. 

The line elements used are superposition of models 
for axial elongation, torsional rotation, shearing and bending. 
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Where stiffeners provide the resistance to bending, the 
classical bending element is used in the analysis. When 
using a facet, the shear bending element is used to provide 
the necessary stiffness from the supporting frames £is this 
provides deformations which are consistent with those of 
the facet. 

As in the facet element, assumed displacements 
are used in the development of the stiffness and loading 
coefficients. The stress coefficients are obtained directly 
from the stiffness matrix. 

In the line element, the local x axis is considered 
to go down the centroid of the beam. The program automati- 
cally makes this transformation if the line element grid 
points are given in the common coordinates. 

Due to the versatility of the program, the input 
necessary from the analysis for the SAMIS porgram is 
extensive. The first section of input for the program is 
the program control cards called "pusedo instructions." 

These cards control tape assignments, matrice maraing and 
ma trice manipulation. Card format is divided into alpha- 
numeric names of the matrices, tape assignment instructions, 
instruction sequence number and program control information. 

The card format is broken up into ten fields. 

The "0" field is the instruction sequence number. The "A", 
"B", and "C" fields contain the tape storage assignments, 
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the fields numbered 1, 2, and 3 contain matrix names, the 
"code" field contains the subroutine or subprogram name and 
the "E" field contains the control information for the 
operations to be performed. Card format is as follows: 



F8.1 


3X, 15 


2X , A3 , 13 


3X , 15 


2X, A3 , 13 


4X, A4 


3X, 15 


2X,A3 ,13 


2X, 16 


0 


A 


1 


B 


2 


CODE 


C 


3 


E 



Psuedo instructions are: 
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TABLE 1 



(Ref. 16, p. 44) 
OPERATION PSEUDO INSTRUCTIONS 



Ins tructions 

ADDS 

BILD 

CHIN 

CHOL 

CONT 

CODE 

COLS 

DECO 

FILL 

FLIP 

INKS 

MULT 

READ 

ROWS 

ROOT 

SAVE 

SORT 

SUBS 

WASH 



Interpretation 

Form C3 = A1 + B2 

Construct small deflection stiffness, stress, 
loading, and/or mass matrices as Al, Bl, Cl, 
and Fl, respectively. 

Form B2 such that B2 B2 T = Al and C3 = B2 _1 
where Al is symmetric and positive definite. 

Form C3 = Al ^ B2 using Choleski decomposition. 

Continuation Card (see pg 43) 

Transform Al to coded format as C3 

Put the Al matrix in column listing and call 
it C3. 

Transform Al to precoded format as C3 . 

Read Al, B2, and/or C3 into core. 

Form C3 = A1 T 

Print matrices Al, B2, and/or C3 
Form C3 = Al B2 

Read matrices Al , B2, and/or C3 from cards. 

Put the Al matrix in row listing and call it C3. 

Find latent roots and vectors of Al, a 
symmetric matrix. Let vectors be B2 , roots C3 . 

Write Al, B2, and/or C3 on tape. 

Sort a matrix Al by row or column as C3 . 

Form C3 = Al - B2 

WASH Al elements from B2 to produce C3 . 
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Psuedo instructions called logic instructions 
provide the capability of loop operations and are also 
required for transfer of control to the success exit. When 
logic instructions are used, the "E" field indicates the 
number of times the logic instruction is to be carried out. 
Logic instructions are: 
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TABLE 2 

(Ref. 16, p. 41) 





LOGIC PSEUDO INSTRUCTIONS 


Instructions 


Interpretation 


PREP 


Prepare for multiple execution of the 
following instructions. Execute the 
instructions between the PREP and the 
next BACK instruction the number of 
times specified in the "E" field. 


VARY 


Vary matrix or tape numbers in the next 
instruction by augmenting corresponding 
field data by the specified integers 
after one pass. 


BACK 


Back up and repeat instructions after 
PREP. 


ERRS 


Disrupt errors can be corrected as 
follows . 


SKIP 


Skip the next "n" pseudo instructions 
where "n" is specified in the !! E” field. 
(Skip cannot be included between a PREP 
and BACK instruction) . 


STOP 


Stop this case and go to HALT. 


PAWS 


Pause in the calculations. Operator can 
restart at any time. (Not operative on 

IBM 7094-7040~ DCS) . 


HALT 


Halt and indicate a successful exit. 
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Control cards needed to direct computations for the most 
general case of static and pseudostatic loading and normal 
modes is a single analysis would be: 
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FORTRAN STATEMENT 



0.0 



1.0 


9001 


WTR1 


9002 


2.0 


9004 


KER2 


10001 


3.0 


9003 


KERl 




4.0 


12001 


KTRl 


9002 


4.1 






11001 


5.0 








5.1 


11001 


DURl 




5.2 


9001 


RBC999 




5.3 








6.0 


9001 


WTR1 




6.1 






11002 


7.0 


11002 


WUR1 




8.0 


11001 


KURl 


12003 


9.0 


11002 


WUR1 


9002 


10.0 


12003 


DCC1 




11.0 


11003 


DDC1 


11004 


12.0 


11004 


VDCl 




13.0 


11002 


WUR1 


9002 


14.0 


12002 


ATC1 


12003 


15.0 


9001 


WTR1 


11005 





CHEX 






ARC1 


READ 


9003 


KERl 


SCR2 


DIED 








ADDS 


12001 


KTRl 


ARC1 


CIIOL 


12002 


ATC1 


KURl 


CONT 








ERRS 








CHOL 


9001 


RBC9 




INKS 








STOP 








CHOL 






WUR1 


CONT 








FLIP 


12003 


WLC1 


WLC1 


CHOL 


9002 


TECl 


TECl 


MULT 


12003 


DCC1 




DECO 


11003 


DDC1 


VDCl 


ROOT 








CODE 


9002 


VCC1 


VCC1 


CHOL 


12003 


DWC1 


DWC1 


ADDS 


11005 


AUC2 


AUC2 


MULT 


9002 


DFC1 



Fig. 1 - Facsimile of Operational Control Cards 
(including normal modes) 

(Ref. 19, p. 9 & 10) 



1 

-1 

-^05 

P 

510 N 4 
920 

N 4 

900 



890 
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FORTRAN STATEMENT 



16.0 




FOCI 






READ 






-1 


17.0 




FOCI 


9002 


DFCl 


ADDS 


12001 


FOC2 




18.0 


11001 


KUR1 


12001 


FOC2 


CHOL 


9002 


DTCl 


910 


18.1 


12002 


DDC1 






CONT 








19.0 


12002 


DDC1 






COLS 


11006 


GCCl 




20.0 


12001 


FOC2 






FLIP 




FOR2 




21.0 




FOR2 


12002 


DDCl 


MULT 


11007 


GCC2 




22.0 


9001 


VJTRl 


9002 


DTC1 


MULT 


12003 


FOC3 




23.0 


12003 


FOC3 






FLIP 




FOR3 




24.0 




FOR3 


9002 


DTC1 


MULT 


11008 


GCC3 




25.0 


110 0 6 


GCCl 






DECO 


12002 


DDCl 


3 


26.0 


12002 


DDCl 






INKS 






305 


27.0 


9002 


DTC1 






FILL 








28.0 


10001 


SCR2 




DTCl 


MULT 


11009 


RCC2 


N 


29.0 


11009 


RCC2 






DECO 


12001 


RDC2 


N 


30.0 


12001 


RDC2 




' 


INKS 






N06 



31.0 HALT 



NOTE: Replace N by the number of structural elements and P 
by N + 1 = (16 ) . 

Replace by the number of cards of codes associated 
with prescribed displacements (.12) . 
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Statement number 9, for example, means take 
matrix WURl on tape 11, position two and multiply it by 
matrix TEC1 on tape 9, position two and place the resulting 
matrix on tape 12, position 3 and call it DCC1. Instruc- 
tions are performed in the sequence that they are placed in 
the input. Any similar type program can use the same 
sequence. Comment cards may be included in the pseudo 
instructions . 

The next section of input defines the model weight 
at each grid point which is desired, in addition to that 
automatically modeled by the SAMIS program from given data. 
The first digits of the row and column numbers correspond to 
the grid point number, and the last digit gives the corres- 
ponding degree of freedom the weight is associated with, i.e. 
1, 2, and 3 the x, y, and z direction respectively; 4, 5, and 
6 the off center weights about the x, y, and z axis. This 
matrix is used to introduce concentrated loadings. Card 
format is as follows: 

Card one 



A3, 13 


16 


16,16 


16 


16 


16 


16 


16 


Matrix 
name & 
no . 


No . of 
cards of 
input 


Matrix 
^ size 


Row(-l) or 
column (1) 
listed 


Precoded (1) or 
coded (0) 



Card two 



16 


16 


16 


16 


16 



Column codes followed by row codes 
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Card three 



E12 . 0 


E12 .0 


E12.0 


Matr.ic 


values row 


or column listed 



For coded format 
Card two 



16 


16 


El 2.0 


16 


16 


El 2 . 0 


Row or column or 
column and row 
listed 


Value or 

matrix 

element 









The next matrix data input are the nonzero 
accelerations associated with the restraints that are not 
modeled by the SAMIS program. Corresponding to the assigned 
loading conditions, gravity (g) accelerations are supplied 
in the x, y, and z direction, and (g)/in accelerations 
forces are supplied about the x, y, and z axis. This matrix 
is utilized to make equilibrium checks. Card format is the 
same as the weight matrix. 

The next section of input gives the nonzero 
elements of the stiffness matrix that are not part of the 
elements listed. This matrix gives the holonomic boundary 
conditions imposed. Card format is the same as for the 
weight matrix. 

Two cards giving the modulus of elasticity and the 
modulus of rigidity and a card of zeroes complete this 
section of input. 
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Tabulated grid point coordinates are then given 
in ascending order. If the coordinates are given in the 
common coordinate system, the SAMIS program makes the 
necessary transformations to local coordinates for genera- 
tion of the various matrices. One grid point is given on 
each card with the following format: 



11 


13 


3X,E7 . 0 , E7 . 0 , E7 . 0 


Sard No. 


Grid point no. 


Grid point coordinates 



Following this section of input is the element 
data cards. This input describes each element as to its 
type, position, size, geometry and the material type it is. 
The card format for each of the elements is as follows: 
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^J 1 



I — I 



CO 






W 



A 

CQ 

cXl i — I 

Eh 



CD 
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TABLE 3 CCont'd) 

(Ref. 16, p. 114) 
ELEMENT DATA FORMAT FOR FACET 



AiiLNaai 

iNawana 

aaawnN 

iNawaaa 

•on QHVO 



Material 

Identity 












Coord . 
Identity 

i 




Coord . 
Identity 


co 


















CO 
















-P 


CD 
















P 










-P 




4-> 




•H 


M 








P 




P 




O 


o 








*P co 




•rH CO 




Pi CO 


•H 








O N 




O csj 




CN) 










Pi 




Pi 




d 


















•rH 










nd 




d 




P 










•H 








■O 










H 




P 














o 




O 




o 


















•rH 




0) 






1 — 1 




• 




-p 


•P 


0) o 






rd 




CO 




CO 


Xi 


• tn !-i td 






O 




d 




rd 


tn 


Oi C, (U MH 






o 




p 




• — 1 


■H 


g «J ^ 






PI CO 




CO CO 




W co 


0) 


0) -Cl o o 






X 








>H 


& 


EH o Ll cn 






‘d 




d 




d 










P 




p 




P 










•H 




•H 




•H 










►P 








d 










H 




EH 




■ EH 




(D 






« 




• 




• 


<D co 


p 






d 




d 




d 


-P 


p <D 






P CO 




P co 




P CO 


p • 


-P tr> p 0) 






O X 




O X 




O X 


-P O 


<tf P CD O 






O 




O 




o 




P cd Pi rd 






O 




O 




o 


■P 


Q) ,n* OiiH 
















(0 CD 


CliUD P 


















E 3 
















p o 
cn ^ 


CD CO 

Fh 






i i 


i 1 


1 . . 


i 

i 

i 


ell J 

: i 

i 

i i 


i 

1 

1 

1 


1 


i 

i i 

i 


1 

1 

1 


I 

1 

i 


i 

i 


I 

1 

1 


i — 1 


CN 


co 




LO 


CO 


r- 


00 


oa 
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d) 


















rd 


p 










-P 








(1) 


rH p 










P 








P CD 


rd -P 










•H 




-P 




^ O 


•H rd 










O 




P 




P 


P P 






-P rH 




P-l rH 




*H i — 1 




P O 


CD CD 






P IS] 




ES] 




O csi 




rd Cm 


-P Q< 






•H 




d 




Pa 




CD 


rd E 






O 




•H 








44 >H 


S a) 






Pa 




P 




d 




CO 


Eh 










o 




•H 












d 








P 
















CD 




o 




H 








P 




-P 
















o 




P 




o 




MH 


CD 










-P 




•H 




O 


P 






I — 1 




•H 




-P 




P 


P 






rd 




-P 




CO 




rd -P 


if) 






O rH 




CO i — 1 




rd i — I 




r-H P 


CO 






O X 




rO X 




rH X 




O (D 


CD 






PI 




P 




W 




E 


p 










CO 






g 


o 


PQ 






d 








d 


22 


si 








P 




-p 




p 


m 










n 




CO 




o 


CQ 










O 




p 




o 




1 — 1 








CD 




•H 




CD 


Pi 


rd 








CO 




Cm 




CO 


— O 


p 


-p 
















in p_, 


I o 




s 


g 


• 


g 


♦ 


g 


• 


i — I 


Cf>H 


tn 






d 


2! 


d 


22 


d 


rH Eh 


(0-P rd 


•H 


w 


H 


P rH 


w 


P rH 


w 


p i — i 




GO Cl) 


CD 


PQ 


CQ 


o x 


CQ 


O X 


CQ 


O X 


♦ § 


^(1) ^ 


s 






o 




o 




o 


Oj Pi 


U02 rtj 








o 




o 




CJ 


w o 




















PI *• Cm 




















PQ co 


ro 


















^ i — I i<] 




















Eh Eh 


O • 


















• <i 


•i-l O 




Pi 


UJ 




g 




g 


-p 


MM Q 


-P g 




Eh 


Eh 




Eh 




Eh 


p 


CD 


(0 




H 


H 


-P rH 


H 




H 


•H i — 1 


Pi Eh 


rd CD 




j<: 


i-4- 


P CS] 


r^» 

l-4» 




r->- 

1-4- 


O CS] 


w g 


rH d 








•H 








Pi 


CQ 


w o 








O 










g 


g 








Pa 








d 


w 


















•H 


PI 










d 








P 


W 


CM 


XiCN 






*H 




CD CM 




o 






-P 






P 




-P 








U • 


■H • 


Q 


Q 


O 


Q 


P • 


Q 


o 




•H O 


P O 


W 


H 




W 


-P O 


W 


•H 




-P g 


P g 


CO 


CO 


i — 1 i — 1 


0) 


*H g 


CO 


-P rH 




CO 


k P 


D 


D 


rd X 


D 


-P 


D 


(0 X 




rd CD 


-P CD 






O 




CO CD 




rd 




i — 1 d 


P d 






O 




Pi r cd 




i — 1 




« o 


o o 






PI 




P O 




W 




g 


o is 










co g 
























i \ 












-M 

CO 








-M 

CO 




r-H 


X rH 






p 




CD rH 




p 






-P 






•H 




P 




•H 




O * 


*H * 










P • 




Cm 




•H O 


P o 










-P o 








-P g 


P g 


Eh 


Eh 


• 


Eh 


‘H X 


Eh 


• 




CO 


*H 


O 


O 


d rH 


O 


-P 


O 


d rH 




rd CD 


-P <D 


g 


g 


P X 


g 


CO CD 


g 


P X 




rH n a 


P d 






o 




Pi d 




O 




w o 


O O 






o 




P o 




o 




g 


O g 






CJ 




co g 




CJ 


xniNaai 




















iNawana 




















H3HWHN 








: I 












iNawaia 




















•ON QHVD 


rH 


CM 


CO 




in 


CD 




co 


o> 
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u 
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B 
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1 /\o_UIb u 

(Ref. 16, pgs . 124, 125, 126 & 127) 
FACET ELEMENT DATA 



Card Columns 




1 


1 


11 


1 


1 


2-4 


13 


K 


1 


5-7 


IX, 12 


31 


1 


8-14 


E7 .0 


GA 




15-21 


E7 .0 


GB 




22-28 


E7 .0 


GC 


1 


29-35 


E7 . 0 


— 


1 


36-42 


E7 .0 


SA 




43-49 




SB 




50-56 




SC 


1 


57-63 


E7 .0 


M 


1 


64-70 


E7 .0 


T 



Indicates card number = 1. 

Element number. 0 < K <_ 9 99. 

Specifies the Facet element 
assumptions. The first digit, 
3, indicates three grid points 
are required. 

First, second, and third 
elastic grid point numbers for 
the element. If negative, 
solution displacements are 
produced in the local coordin- 
ate system at the grid point. 

Data in this field will be 
ignored . 

First, second, and third sub- 
stitute grid point numbers . 

If the grid point numbers are 
negative, solution displace- 
ments are produced in the local 
coordinate system at the grid 
point . 

Facet mass ^ 

M <_ 0 , total mass:#sec /in 
M > 0, mass per unit area: 
ftsec^/in-^ 

Facet thickness in inches. 



1 71-72 A2 N Name of the structural material. 

The first two characters of the 
material name must match the 
first two characters of the 
material name in the material 
table . 

Numbers given in this column are to be taken literally. The 
user must substitute appropriate numbers for letters given. 
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TABLE 5 (Cont'd) 



FACET ELEMENT DATA 



* 

Card Columns Format Item Interpretation 



2 


1 


11 


2 


2 


2-4 


13 


K 


2 


5-7 


IX, 12 


31 


2 


8-14 


E7 .0 


CA 




15-21 


E7 .0 


CB 




22-28 


E7 .0 


CC 


2 


29-35 


E7 .0 


— 


2 


36-42 


E7 .0 


P 


2 


43-49 


E7 .0 


T 

m 


2 


50-56 


E7 .0 


T u 


2 


57-63 


E7.0 


T 1 


2 


64-70 


E7 .0 


— 


2 


71-72 


A2 




3 


Omit 






4 


Omit 







Indicates card number = 2. 

Element number. 0 <_ K 999. 

(may be omitted) 

Same as Card 1, columns 5-7. 

(may be omitted) 

Continuity boundary condi- 
tions at the first, second, 
and third grid points. If both 
substitute points are given, 
these apply to the substitute 
points . If only elastic appear, 
they apply to the elastic . 

Ignored 

2 

Normal pressure: pounds/inch 
Positive in the plus z 
direction . 

Temperature (Degrees Rankine) 
of the material (used to de- 
fine elastic constants) . 

Upper surface temperature 
change (Degrees Rankine) from 
zero stress temperature. 

Lower surface temperature 
change (Degrees Rankine) from 
zero stress temperatuer . 

Ignored 

Ignored 



Numbers given in this column are to be taken literally. The 
user must substitute appropriate numbers for letters given. 
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TABLE 5 (Cont'd) 
FACET ELEMENT DATA 



* 

Card Columns Format Item Interpretation 



5 


1 


11 


5 


5 


2-4 


13 


K 


5 


5-7 


IX, 12 


31 


5 


8-70 


9E7 .0 





5 


71-72 


A2 


— 


6 


Omit 






7 


1 


11 


7 


7 


2-4 


13 


K 


7 


5-7 


IX, 12 


31 


7 


8-19 


9E7 .0 




7 


71-72 


A 2 


C 



3 Omit 



Indicates card number = 5. 

Element number. 0 £ K £ 999. 
(may be omitted) 

Same as Card 1, columns 5-7. 
(may be omitted) 

Coordinates (inches) in the 
overall system of the origin of 
the local coordinate system 
(xi /Y i /Z i) ; a point on the 
local x axis {X2,Y2' 7 ‘2^ > an< ^ a 
point in the x-y plane (* 3 ,y 3 , 
Z3), noncollinear with the 
first two points are selected 
to define desired direction of 
the local z axis. 

Ignored 



Indicates card number = 7 . 

Element number. 0 <_ K £ 999. 
(may be omitted) 

Same as Card 1, columns 5-7. 
(may be omitted) 

Coordinates (inches) of the 
substitute grid points corres- 
ponding to the first, second, 
and third elastic grid points. 

Coordinate identification. If 
C = L, substitute grid point 
coordinates are in the local 
system. If C / L, coordin- 
ates are overall . 



Numbers given in this column are to be taken literally . The 
user must substitute appropriate numbers for letters given. 
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TABLE 5 (Cont'd) 



FACET ELEMENT DATA 



•k 

Card Columns Format Item Interpretation 



9 


1 


11 


9 


9 


2-4 


13 


K 


9 


5-7 


IX, 12 


31 


9 


3-70 


9E7 . 0 





9 71-72 A2 C 



Card number = 9 . 

Element number. 0 <_ K £ 999. 

Same as Card 1, columns 5-7. 

(may be omitted) 

Coordinates (inches) for the 
first (x^ ,y]_ / Z]_) , second (X 2 , 

y 2 / z 2 ) and third (*3^3,23) 

elastic grid points of the 
Facet in the overall or local 
system . 

Coordinate identification. If 
C = L elastic grid point coor- 
dinates are in the local system. 
If C ?! L, coordinates are over- 
all . 



Numbers given in this column are to be taken literally. The 
user must substitute appropriate numbers for letters given. 
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rn 7\ "D T TP C 

O. /II-JJ-JJLJ KJ 



(Ref. 16, pgs . 128, 129, 130, & 131) 
2 LINE ELEMENT DATA 




1 


1 


11 


1 


1 


2-4 


13 


K 


1 


5-7 


IX, 12 


2A 


1 


8-14 


E7 . 0 


GA 




15-21 


E7 .0 


GB 




22-28 


E7 .0 


GC 


1 


29-35 


E7.0 


A 

X 


1 


36-42 


E7 .0 


J x 


1 


43-49 


E7 .0 


A 

y 


1 


50-56 


E7 .0 


i 

z 



Indicates card number = 1. 

Element number. 0 < I< £ 999. 

Specifies the line element equa- 
tions . Rod and tube equations 
and (A = 1) elementary beam or 
(A = 2) shear beam displace- 
ment assumptions . 

First, second, and third grid 
point numbers for the element. 
The grid number of the third 
grid point may be omitted. Tf 
numbers are negative, solution 
displacements are in the local 
coordinate system at the grid 
point . 

Cross-sectional area (axial) 
in the y-z plane: inches . 

Torsional rigidity against 
twist about the x axis: 
inches ^ . 

Effective area deforming in 
shear in the x-y plane due to 
a y force: inches^. 

Moment of inertia resisting a 
moment about the z axis: 
inches^ . 



'k 

Numbers given in this column are to be taken literally. The 
user must substitute appropriate numbers for letters. 
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TABLE 6 (ConL’d) 



LINE ELEMENT DATA 

Card Columns Format I tern* Interpretation 



1 


57-63 


E7 .0 


A 

z 


Effective shear area for a 
shear force in the z direc- 
tion: inches^. 


1 


64-70 


E7 .0 


I 

Y 


Moment of inertia resisting a 
twist about the y axis: 

inches^ . 


1 


71-72 


A2 


N 


Name of the structural mater- 
ial. The first two characters 
of the name must match the 
first two characters of the 
material table material name. 


2 


1 


11 


2 


Card number = 2 . 


2 


2-4 


13 


K 


Element number. 0 < K < 999. 


2 


5-7 


IX, 12 


2A 


Same as Card 1, columns 5-7. 
(may be omitted) 


2 


8-21 


2E7.0 


CA 

CB 


Continuity boundary conditions 
at the first and second grid 
points. If substitute grid 
points are given, these apply 
to the substitute points. If 
only elastic appear, they 
apply to the elastic. 


2 


22-28 


E7 .0 


--- 


Ignored 


2 


29-35 


E7 . 0 


M 


Line element mass. If M < 0, i 
the total mass (Osec^in) . If 
M > 0 , M is the mass per unit 
length (ttsec^/in^) . 


2 


36-42 


E7 .0 


P 


Normal pressure: pounds/inch 
(positive in the plus z 
direction) . 


2 


43-49 


E7.0 


T 

m 


Temperature (degrees Rankine) 
of the material (used to de- 
fine elastic constants). 


2 


50-56 


E7.0 


T 

o 


Mean temperature change (de- 
grees Rankine per unit of 
cross-sectional area - A x ) of 
the element from the zero 
stress temperature. 



Numbers given in this column are to be taken literally. The 
user must substitute appropriate numbers for letters. 



TABLE 6 (Cont'd) 



LIME ELEMENT DATA 



Card 


Columns 


Format 


Item* 


Interpretation 


2 


57-63 


E7 .0 


T 

Y 


Temperature gradient (degrees 

Rankine per unit cross- 

sectional moment of inertia 

I ) in the z direction, 
y 


2 


64-70 


E7 .0 


T 

z 


JL 

Temperature gradient (degrees 

T9 a ri i ra nor n m' +- . 



Rankine per unit cross- 
sectional moment of inertia 
I ) in the y direction. 



2 


71-72 


A 2 


— 


Ignored 


3 


Omit 








4 


Omit 








5 


1 


11 


5 


Indicates card number = 5. 


5 


2-4 


13 


K 


Element number. 0 < K < 999. 
(may be omitted) 


5 


5-7 


IX, 12 


2A 


Same as Card 1, columns 5-7. 
(may be omitted) 


5 


8-70 


9E7 ,0 




Coordinates (inches) in the 
overall system of the origin 
of the local coordinate system, 
(xi,yi,zi) ; a point on the 
local x axis ( X 2'Y2' Z 2) ; anc * a 
point in the x-y plane (^ 3 ^ 3 , 
Z 3 ), noncollinear with the 
first two points and located 
to define desired plus z 
direction . 


5 


71-72 


A 2 


— 


Ignored 


6 


Omit 








7 


1 


11 


7 


Card number = 7 . 


7 


2-4 


13 


K 


Element number. 0 < K < 999. 


7 


5-7 


IX, 12 


2A 


Same as Card 1, columns 507. 
(may be omitted) 



Numbers given in this column are to be taken literally. The 
user must substitute appropriate numbers for letters. 



TABLE 6 (Cont'd) 



LINE ELEMENT DATA 



Card Columns Format Item * 

7 3-21 2E7.0 SA 

SB 



7 


22-28 


E7 . 0 


-- 


7 


29-70 


6E7.0 


— 



7 


71-7 2 


A 2 


C 


3 


Omit 






9 


1 


19 


9 


9 


2-4 


13 


K 


9 


5-7 


IX, 12 


2A 


9 


8-70 


9E7 .0 






9 71-72 A2 C 



Interpretation 

First and second substitute 
grid point numbers . If the 
grid point numbers are nega- 
tive, solution displacements 
are produced in the local 
coordinate system at the 
grid point. 

Ignored 

Coordinates (inches) of the 
substitute grid points corres- 
ponding to the first and 
second elastic grid points: 
all in the local or overall 
system . 

Coordinate identity. If C = L, 
substitute grid point coordin- 
ates are in the local system. 
If C / L, coordinates are 
overall . 



Card number = 9 . 

Element number. 0 <_ K £ 99 9. 
(may be omitted) 

Same as Card 1, columns 5-7. 
(may be omitted) 

Coordinates (inches) for the 
first (x 1 ,y 1 , z^) , and second 
(x 2 ,y 2 r z^) elastic grid points 
of the line element and the 
grid point (x 2 ,y 3 ,z 2 ) defining 
the principal plane of the 
cross-section: all in the over- 
all or local systems. 

Coordinate identification. If 
C = L, elastic grid point 
coordinates are in the local 
system. If C / L, coordinates 
are in the overall system. 



Numbers given in this column are to be taken literally. The 
user must substitute appropriate numbers for letters. 
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For each material used in the structure two cards 
of material data arc required. These cards give the material 
identification, the temperature associated v/ith the material 
properties, the coefficient of material expansion and 
material stiffness coefficients. The material stiffness 
coefficients are the coefficients of the stress-strain equa- 
tions in accordance with the following equations: 



C *^1 

cr 

XX 




D n 














M 


a yy 




°21 


°22 












\y 


j 


T ~~~ 


°31 


°32 


°33 










1 XY r 


\ a zz 




°41 


°42 


°43 


°44 








\ < zz 


Icr 

1 xz 












°55 






\z 


v. yz . 












°6 5 


D 

Oh 

Oh 

1 


1 V 


VZ 

J J 



The above equations hold true for a monotropic material. If 
= D 43 = Dgj = 0 an orthotropic material is des- 
cribed. An isotropic material is described if in addition 

D 33 = °55 = D 66 = D 11 = °41 = °42 and °21 = + ^ U ll^ * 

Then D cc = E/(2(l + v)) and D.. = E(1 - v) / ( (1 + v)(l - 2v) ) 
5 5 11 

Card format is as follows: 
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MATERIAL TABLES INPUT DATA 



o 



o 























o 








CO 


CO 










O MO 




o 




0 




CD 


CD 










• + 


O 


o 




> 1 




•H 


CD 


i — 1 








CO p 




o 


G 


0 cn E 




-P 


P 


CO 








W • 




o 


0 


,g Td p o g 




P 


tr> 


a 


P 










o 


0 


(D 0 £ G 




CD 


CD 
















p 


-P Td XI -P 




04 


Td 














o 




0 G g G 




0 


1 


CM 


co 










o 


0 


G 0 g tn o 




p 


rG 


CM 












o 


> 


E E £ £ -H 




04 


u 


a 


P 






o 


MO 


o 


0 


E -H -P 






G 










• o 


+ 


o 


XJ 


•H O G Td 0 




1 — 1 


•H 




•w 






oo 


P 


o 




0 O 0 0 u 




fd 


\ 


1 — 1 


CM 






W 


• 


o 




•H 0 -H 0 -H 




•H 


CO 


CM 








^ — 




o 


0 


P P -P rH MH 




P 


CD 


a 


P 












1 — 1 


<D (d -H 




CD 


XJ 














o 


0 


■P CO *H (D -P 




-P 


o 














o 


•H 


td -H O rC g 




fd 


G 


• — i 


i — i 










o 


p 


E 0 -P 0 




E 


•H 


i — i 








o 




o 


0 


-P co Td 








a 


P 






• o 


O 


o 


-p 


rG h cn >»-h 




G 


• • 










CO 




o 


0 


O < rH 




CD 


G 




« 






W 




o 


E 


(d • a p 




> 


O 


CO 


co 






— " 




o 




W P E O 0) 




■H 


co 


-p 


-P 












1 — 1 


O 3 -P 






•H 


G 


G 










o 


1 — 1 


•• X) G * o 






G 


CD 


CD 










o 


0 


G E *rH d) fd 




CD 


td 


■H 


*H 






^ MO 




o 




0 G E rH P 




X! 


04 


O 


O 






o H- 


MO 


o 


p 


•h d d X) d 




-P 


X 


•H 


•H ♦ 






• W 


+ 


o 


0 


■P i — 1 # H 






CD 


MH 


MH G 






co • 


P 


o 


MH 


rd G <1 cn o 




P 




MH 


MH -H 


f=H 




W vo 


• 


o 




o 0 Ml 




0 


i — 1 


CD 


0 • 


b 




^ rH 


H' 


o 




•H -H d O X 




MH 


fd 


0 


0 O 1 


p 










0 


4-> G CVH 


# 




E 


O • 


O co 


p 








o 


Td 


•h id (d co 


-P 


CD 


p 


G 


\ 


H 


E 






o 


p 


■P o co 


G 


P 


CD 


CO -H 


CO CO 




G 




MO 


o 


0 


G *H P fd 0 


rd 


G 


gj 


CO * 


co G 


w 


G 


O MO 


+ 


o 


o 


a) mh <; rC 


O 


-P 


-p 


CD tJ 1 


0 rH 


p 


•H 


• + 


P 


o 




Td *H CO p -P 


-H 


fd 




G CO 


G 


p 


E 


CO W 


* 


o 


0 


•h -P rd 


MH 


P 


MH 


mh \ 


MH - 


< 


G 


H • 


MO 


o 


1 — 1 


G Td MH 4-1 


*H 


CD 


0 


MH CO 


MH MO 


Eh 


rH 


/ CO 


i — 1 


o 


x> 


rH CD P 0 0 


G 


04 




•H ,.Q 


•H MO 




< 








0 


fd td co 


tP 


E 


-P 


-P rH 


-P P 


P 








o 


-p 


*H *(—} Tj G CO 


•n 


<L> 


rH 


CO 


CO 




o 






o 




P G -H p 


CO 


-P 


CD 


■K 






•H 


^ MO 




o 


1 — 1 


0 0 0 CD 






•H 


1 — 1 CO 


i — i in 


Pd 


CP 


O + 




o 


0 


D d Dd D 


0 




O 


rd co 


0 MO 


w 


o 


• W 


O 


o 


*H 


fd tT 1 CO CD O 


p 


G 


•H G 


■h a 


•H p 


Eh 


p 


CO • 




o 


p 


E -h co fd 


fd 


•H 


mh *h 


p 


P 




-p 


W vo 




o 


0 


G -P G P 




Pd 


mh m 


CD - 


0 


X 


0 


^ i — i 




o 


-p 


(D g fd <d 


P 


G 


CD G 


•P CM 


-P ID 




0 








0 


Xd GOG 


(D 


<d 


0 (d 


(d co 


0 m 


W 


M 






o 


E 


EH RJ +> 43 OX! 


Pd 


U Pd 


S Q 


S P 


P 


— 






o 
















P 




1 




o 


MH 














s 




o w 


MO 


o 


0 


















• in 


+ 


o 
















CO 




CO CM 


P 


o 


U) 


















W rH 


• 


o 


p 


















— " • 


CO 


o 


♦H 
























0 






















o 


p 


-P 


C 








, — . 


<-> 










o 




fd 






o 


o 


o 


o 






— S CT) 




o 


1 — 1 


E 






• 


• 


• 


# 






o H- 


MO 


o 


1 — 1 


p 


< 




CO 


00 


CO 


CO 






• W 


+ 


o 


0 


o 


»«. 




W 


w 


W 


p 






CO CO 


W 


o 






X 




— 


v — 


MO 








P cn 


• 


o 


-p 


CM 








— " 


— ' 






^ • 


CO 


o 


0 
























XJ 
























-p 










































^ vo 




o 


0 


















- Eh 




o 


-p 


















cm ^ 




o 


0 


Td 


















fG i — l 




o 


o 


1 — 1 










CTi 


CO 






o 




o 


*H 


CD 


rH 




CM 


co 


1 


1 






X CM 




o 


Td 


•H 










h 1 


CM 






CM 




o 


G 


P 






















o 


*H 
























0 






















* 


-P 


















Eh rH 


CV] 


CO 




















£5 






Td 


Td 


















X P 


p 


p 


0 


p 


rH 




rH 


rH 


i — i 


CM 






K K K 


Pd 


U) 


fd 


















O << «< 




D 


CJ 


















IL U 


u 


o 


* 
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The material tables must be followed by one card of zeroes 
to signal the end of the material tables. 

The next input is the nonzero boundary conditions 
applied to displacements or loads . Displacements are given 
in inches for the restrained degrees of freedom in the x, y, 
and z directions and in radians for the restrained degrees 
of freedom about the x, y, and z axis. Loads applied in the 
x, y, and z directions are given in pounds and those applied 
about the x, y, and z axis are given in inch-pounds. The 
loads are those associated with the unrestrained degrees of 
freedom. Card format is the same as that for the weight 
matrix . 



The last section of input is the title information 
cards. There are a total of eleven cords defining the dis- 
placements, generalized stiffness and generalized weight 
matrices. This is written as follows: 



DDC 1,2,3 = DISPLACEMENT, GENERALIZED STIFFNESS AMD GEN- 
ERAL WEIGHT MATRICES 

ROW AND COLUMN CORRESPOND TO MODE, LOAD, OR REACTION 
DISPLACEMENT OR EXCEPT DDC L ROW LEADING DIGITS CORRESPOND 
TO GRID POINT AND FINAL DIGIT TO THE DEGREE OF FREEDOM AT 
THE GRID POINT 

REACTION NOT DISPLACEMENT AT RESTRAINT 

INTERNAL FORCE FOR ELEMENT -X- (RDC X) 

ROW LEADING DIGITS CORRESPONS TO GRID POINT, FINAL DIGIT 

1 AXIAL FORCE 4 TORQUE 

8 SHEAR ALONG X2 5 MOMENT ABOUT X2 

7 SHEAR ALONG X3 6 MOMENT ABOUT X3 

COL CORRESPOND TO MODE, LOAD, OR REACT. DISP. OR ACC. 
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The complete computer program is available in the 



Naval Architecture and Marine Engineering Library. 
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MODELING OF TIIE STRUCTURE 



The structure analysis of a system as complicated 
as an icebreaker bow requires considerable preparation in 
order to use the SAMIS program. The first requirement is to 
develop a three-dimensional picture of the structure and 
orient it in relation to the common coordinate system. In 
order to clearly indicate the necessary details and label 
the elements and grid points, the scale must be quite large. 
In Appendix A, the three-dimensional drawing of the ice- 
breaker VJESTWIND is in a scale of 1/2" = 1'. Once the 
three-dimensional table has been drawn, a material table 
must be set up to give the necessary material parameters 
required for the input data. Appendix B gives the material 
table for the icebreaker VJESTWIND. 

Since a ship's structure is symmetrical about 
the centerline, only one half of the ship need be drawn and 
analyzed. The boundary conditions for the nodes on the 
centerline can be completely described by displacements . 
Unless actual whipping of the ship occurs, all nodes can 
be fixed in the athwartship direction. Movement in the fore 
and aft direction will be sufficiently modeled by the facets 
edge, and movement in the vertical direction can be modeled 
as very stiff spring. 
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The boundary conditions of the hull are simply the 
average pressure loads for the depth the facet is beneath 
the waterline. At the bow, shell pressure loads for the ice 
contact area can be selected as desired. The total area in 
contact with the ice is that which will give the load res- 
ponse determined by Dr. White (Ref. 5). Shear and bending 
moment diagrams developed from the loading conditions imposed 
can be used to give the boundary forces where the structure 
is cut off. These can be applied as concentrated loads in 
the matrix data as was described as Matrix WTRl in the 
previous section. This will necessitate some hand calcula- 
tions to determine the section modulus and shear area. 
Furthermore, it would be necessary to determine the 
deacceleration at the time of maximum bow load and apply 
that to the mass of the remainder of the hull to get a 
deacceleration force. A check is needed to determine if 
these total forces satisfy equilibrium for the entire 
system being analyzed. A simpler method would be to assign 
large stiffnesses to springs connected to fixed points and 
supply these stiffness as matrix data input as was described 
by matrix KERl in the previous section. 

It would be desirable to divide the entire section 
being analyzed into small facets and model each stiffener as 
a beam. However, due to the limit on the number of elements 
the program can handle (999), this type of modeling can only 
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be done in those areas v/here complete and accurate detail 
is desired. The IIT Research Institute Technology Center 
of Chicago, Illinois in IITRI Project J6127 (Ref. 19), 
developed what is called an orthotropic plate model for 
stiffened structures so that large facets may be used for the 
remainder of the structure. Thickness is modeled as satis- 
fying the following equation. 

2 12 I 

t A* + A 

f P 

where 



A 



P 

■f 

I 

d 



bh (the plate area 
the area of the frame 



A 



A c d‘ 

f 



P 



A + A, 



+ I + I, 

P J 



distance between the centroids of the plating 
and frame. 



If = Moment of inertia of frame alone about its own 
centroidal axis. 



I = ~- 

p 12 (1-v ) 

b = distance between frames 
h = plate thickness. 



The modulus of elasticity in the plane of the plate and in 
the direction the stiffeners run is 
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E 



E (A r + A ) 

= f P 

y 2 I p 

bt(l-v 



The modulus of elasticity iix the plane of the plate and 

• * 

perpendicular to the plate is 

I 

E = E 

x I y 

and 



E = vE 
xy x 



The resulting material constants are then 



D, , = E + E /E 
11 y xy' y 

D 0 , = E + E 2 /E 

21 xy xy' y 

0 oo = E + E 2 /E 

22 x xy' y 



°31 = 0 



°32 = 0 



°33 G 



°41 = E xy 
°42 E xy 



°43 = 0 



°44 = E y 



D 55 = G 
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0 



D 



65 



D 
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G 



where G 




Utilizing the SAMIS program, version one IITRI analyzed a 
section of the WESTWIND hull from frame seven to thirty-one. 
Results were not compared with stress data accumulated by 
the WESTWIND, but showed a structural weakness at frame 25 
which has been an area of structural failures in the past. 



66 



FUTURE DEVELOPMENTS 



The capabilities of this program provide the ships 
structure analysts with a very valuable tool. Modeling the 
entire shell of a ship as an unstiffened plate, and then 
designing the stiffeners to be most effective from the 
results of the first run can lead to optimal utilization of 
material in ship construction not only for tough skinned 
vessels such as an icebreaker, but also for all other types 
of vessels. At last there is a glimmer of hope that we can 
depart from the long standing method of using something that 
worked before in a new ship so we can have hope of some 
measure of success. 
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APPENDIX B 



MATERIAL TABLES 

(All material is medium black steel unless otherwise indicated) 

Shell Plating 

Main deck to second deck. 

20.4# high strength steel plate. 

Second deck to platform deck. 

66.3# high strength steel plate. 

Platform deck to keel. 

56.1# high strength steel plate. 

Transverse Frames 

Stem, Main deck to 30° bow. 

15" x 15.3# plate. 

30° bow, frame 7 to keel. 

40" x 20.4# plate. 

30° bow, frame 7 to keel. 

20 " x 40.8# plate . 

Cant frame 02 to frame 61, main deck to second deck. 

6" x 4" x 8.25# T 
Second deck to keel. 

15" x 4" x 15.3# T 

Note special framing detail in Appendix C. 
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Transverse Bulkheads 



(All stiffeners on 24" centers) 
Bulkhead 7 
Plating 
7.65# plate. 

Stiffeners 
6" x 4" x 11# T 
Bulkhead 19 
Plating 

Main deck to second deck. 

7.65# plate. 

Second deck to platform deck. 
12,75# plate. 

Stiffeners 

6" x 4" x 3.25# T 

Bulkhead 31 

Main deck to second deck 
7.65# plate . 

Second deck to platform deck. 
10.2# plate. 

Platform deck to keel. 

12.75# plate. 

Stiffeners 

Main deck to second deck. 

6" x 4" x 11# T. 
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Bulkhead 43 



Plating 

Main deck to third deck 
7.65# plate. 

Third deck to keel . 

10.2# plate. 

Stiffeners 

Main deck to second deck 
5" x 2.69" x 4.48# T 
Second deck to shell. 

6" x 4" x 11# T 
Bulkhead 61 
Plating 

Main deck to third deck. 
7.65# plate. 

Third deck to keel 
10.2# plate. 

Stiffeners 

Main deck to second deck 
5" x 2.69" x 4.48# T 
Second deck to shell 
6" x 4" x 11# T 
Longitudinal Bulkheads 



Frame 31 to Frame 61 



Main deck to second deck. 



10.2# plate 

Second deck to shell. 

12.75# plate 
Stiffeners 

Frame 31 to frame 61 on 16" centers 
6" x 4" x 11# T. 

Decks 
Main deck 
Plating 

Stem to frame 23 
10.2# plate. 

Frame 28 to frame 61. 

7.65# plate. 

Transverse Stiffeners 
Stem to frame 61 
6" x 4" x 8.25# T. 

Frame 34, from centerline to 10' on either side of centerline. 
12" x 4” x 16.5# T 
Longitudinal stiffeners 

Frame 7 to frame 19, 1' off centerline at frame 7 and 5' off 
centerline at frame 19 . 

12” x 6 1/2” x 25# I 

Frame 31 to frame 61, 9' off centerline at frame 31 and 19' 
off centerline at frame 61. 
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12" x 6 1/2" x 25# I 
Second deck 
Plating 

Bow to frame 31 
7.65# plate. 

Frame 31 to frame 61 
10.2# plate. 

Deck Stiffeners. All 
6" x 4" x 11# T. 

Transverse beam frame 35 
12" x 4" x 25# T. 

Longitudinal beams 

Bow to frame 7 on centerline. 

12" x 6 1/2" x 25# I 

1' from centerline at frame 7 to 6.0 from centerline at 
frame 27. 

12" x 6 1/2" x 25# I. 

Frame 31 to frame 43 on centerline. 

12" x 6 1/2" x 25# I. 

Frame 43 to frame 61, 6.0' from centerline. 

12 x 6 1/2" x 25# I 
Third deck 
Plating (All) 

10.2# plate. 

Transverse stiffeners 
Bow to frame 61 
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7" x 6 3/4" x 15# T. 

Frame 45 and frame 49 from centerline to 6' on either side 
of centerline. 

12" x 6 3/4" x 30# T 
Longitudinal stiffeners 

1' from centerline at frame 7 to 9 ' from centerline at 
frame 31. 

Frame 31 to frame 43 on centerline. 

Frame 43 to frame 61, 6* from centerline. 

16" x 7" x 36# I 

12.75# brackets, 21" x 10 1/2", used to tie into frame at end 
of each deck stringer. 

Platform Deck 
Plating 

Frame 7 to frame 31 and frame 43 to frame 61. 

10.2# plate. 

Transverse stiffeners (All) 

7" x 6 3/4" x 30# T. 

Frame 44 and frame 48, from centerline to 6 ' on either side 
of centerline. 

12" x 6 3/4" x 30# T. 

Longitudinal stiffeners 

Frame 43 to frame 56 on centerline. 

16" x 7" x 36# I. 
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